Cascades and dissipation ratio in rotating MHD turbulence at low magnetic Prandtl 
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A phenomenology of isotropic magnetohydrodynamic turbulence subject to both rotation and 
applied magnetic field is presented. It is assumed that the triple correlations decay-time is the 
shortest between the eddy turn-over time and the ones associated to the rotating frequency and 
Alfven wave period. For Pm = 1 it leads to four kinds of piecewise spectra, depending on the four 
parameters, injection rate of energy, magnetic diffusivity, rotation rate and applied field. With a 
shell model of MHD turbulence (including rotation and applied magnetic field) , spectra for Pm < 1 
are presented, together with the ratio between magnetic and viscous dissipation. October 18, 2010 
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I. INTRODUCTION 

Magnetohydrodnamic turbulence in natural objects is 
often subject to global rotation or applied magnetic field, 
or both. In the Earth's core the turbulence occurs under 
the fast rotation of the planet and is embedded in the 
dipolar magnetic field produced by dynamo action. Such 
double effect is currently studied in an experiment with 
liquid sodium [l[. Waves of different types have been 
measured that might be attributed to either Alfven or 
Rossby waves or a combination of both. The frequency 
spectra show a scries of bumps, attributed to wave fre- 
quencies, in addition to piecewise slopes. A proper un- 
derstanding of such rotating MHD-turbulence would re- 
quire a non-isotropic formalism. Several ones have been 
developed for fast rotation Phenomenological ap- 

proaches relying on three-wave [5| or four-wave [6j reso- 
nant interactions have been developed for an applied field 
and documented numerically Q- 

In the present paper we come back to the Iroshnikov 
H and Kraichnan g phenomenology for isotropic MHD 
turbulence. They argue that the destruction of phase co- 
herence by Alfven waves traveling in opposite directions 
introduces a new time-scale ta- It might control the en- 
ergy transfer, provided it is shorter than the eddy turn- 
over time-scale tk- Applying the same idea, Zhou [l(| 
suggests that due to global rotation the kinetic energy 
spectrum is affected through phase scrambling, leading 
to a third time-scale associated to the rotation fre- 
quency. The generalization to both global rotation and 
applied magnetic field is therefore straightforward (see 
section |H| , the energy transfers being controlled by the 
shortest time-scale between tk, ta and rn- 

An advantage of assuming isotropy is that it can 
be tested against simulations with shell models. Shell 
models are toy-models that mimic the original Navier- 
Stokes and induction equations projected in Fourier 
space, within shells which are logarithmically spaced. 
There arc only two complex variables per shell, one cor- 
responding to the velocity, the other to the magnetic field 



[111 [l2| . Depending on the model, the energy transfers 
can be considered as local or not [13[ . Such models allow 
for simulations at realistically low viscosity v and mag- 
netic Prandtl number Pm = vjr\ [3] , where r\ is the mag- 
netic diffusivity. The time dependency of the solutions 
is strongly chaotic, eventually leading to intermittency. 
Therefore, though all geometrical details of velocity and 
magnetic fields are lost, shell models give relevant in- 
formations on spectral quantities like energies, helicities, 
energy transfers, etc. In section IIIII we introduce such 
a shell model of rotating MHD turbulence, taking care 
to keep the terms corresponding to rotation and applied 
magnetic field as simple as possible. For Pm < 1 we cal- 
culate the spectra for different values of rotation D, and 
applied field Va- We also calculate the ratio of the joule 
dissipation over the viscous dissipation, which cannot be 
estimated from scaling laws. 



II. PHENOMENOLOGY 

A. Time scales 

Following @ (see also Q3 and [Hj]), we assume that 
for homogeneous isotropic statistically steady turbulence 
the decay of triple correlations, occurring in a time scale 
T3(fc), is responsible for the turbulent spectral transfer 
e from wavenumbers lower than k to higher wavenum- 
bers. This implies T^{k) ~ e. Assuming in addition that 
e depends only on the wave number k and the kinetic en- 
ergy spectral density E(k), a simple dimensional analysis 
leads to 



T 3 (k)E 2 (k)k 4 



(1) 



The kinetic energy spectral density is defined as E(k) = 
k~ 1 u 2 (k) where u{k) is the characteristic velocity of ed- 
dies at scale k. 

In absence of applied magnetic field and rotation, the 
time scale for T^ik) is the eddy turn-over time 



TK (k) = [ku(k)y 



(2) 



2 



leading to the Kolmogorov turbulence energy spectrum 
E(k) ~ £ 2 / 3 fc" 5 / 3 . 

For fully developed MHD turbulence at Pm = 1 the 
same Kolmogorov spectrum is assumed for both kinetic 
and magnetic energy provided that the system is much 
above the onset for dynamo action [k|. In that case E(k) 
denotes either the kinetic or magnetic energy spectral 
density. In presence of an applied magnetic field Bo an 
other possible time scale for r 3 (fc) is the Alfven time scale 

r A (k) = [Wa}~\ (3) 
leading to the Alfven turbulence energy spectrum E(k) ~ 

Finally for rotating turbulence caused by uniform rota- 
tion £1 a third possible time scale for T3 (k) is the rotating 
frequency 

to = n~\ (4) 

leading to the rotating turbulence energy spectrum 
E(k) ~ fiVV' 2 *- 2 . 

The value of Ts{k) is naturally defined by 

r 3 (k) = win {t k (k), r A (k), m} ■ (5) 

It corresponds to the fastest way to transfer energy to 
smaller scales, between non-linear eddy cascade, Alfcn 
waves interactions and phase scrambling due to rotation. 
In addition we define the magnetic dissipation time scale 
by 

r v (k) = (k 2 V )-\ (6) 

The dissipation range corresponds to k > k v with k v 
defined by T 3 (k n ) = T v (k n ). 

Therefore at each scale fc _1 , we have to compare the 
four time scales TK{k),TA(k),TQ and ^(/c) to figure out 
what kind of turbulence occurs. 



B. Spectra for Pm = 1 

At k < 1, tq < mm {T K (k),TA(k),T v (k)} implying 
that T3(k) = tq, unless £7 = 0. This corresponds to 
a rotating turbulence with E(k) = f2 1//2 e 1 / 2 fc -2 . For 
larger k, TK{k),T A (k) and r^(fc) decrease while rn stays 
constant. Therefore, provided that the dissipation is 
not too strong, a first transition occurs at a scale for 
which ro = mm{rK (k), T A (k)}. This scale is, either (i) 
fci = (0 3 /e) 1/2 if £ > VLV%, or (ii) k x = Vl/V A if e < QV 2 . 
This transition leads to, either (i) a Kolmogorov E{k) = 

£ 2/3 fc -5/3 ; Qr ( n ) an Alfy(§n E (Q = yl/2 £ l/2 fc -3/2 tm _ 

bulence. This transition does not occur if the dissipa- 
tion overcomes the Kolmogorov and Alfven turbulence, 
namely if (i) r\ > e/Vl 2 and (ii) ij > V\/Vt. In that case 
the dissipation scale is given by k n = (O//7) 1 / 2 . 

In case (i) provided again that dissipation is not too 
strong, a second transition occurs at k^ = e/V A . 



This transition leads to an Alfven turbulence E(k) = 
y^/ 2 £ i/2^,-3/2 un ^ji t nc dissipation becomes dominant 
for k > k v with k n = Va/t). If V\ < rje the dissipa- 
tion overcomes the Alfven turbulence and the dissipation 
scale is given by k v = e 1 / 4 ?^ 3 / 4 . 

In case (ii) a second transition toward a Kolmogorov 
turbulence is not possible. Indeed, it would occur at k = 
sjV\ which can not be larger than k\ from the condition 
e < VtV A . In that case the Alfven turbulence simply 
extends to the dissipation scale given by k n = Va/ti- 

The four possible types of inertial regimes are sketched 
in Fig. m m which the spectral energy density is plotted 
versus k for Pm = 1. The slopes and characteristic wave 
numbers are indicated. The conditions to get one of these 
four possible inertial regimes are summarized in the plane 
(Va , fi) in Fig. [21 The case without rotation corresponds 
to the abscissa axis. Then two regimes KA or K are pos- 
sible depending whether r\jV\ < V\l e or n °t- The case 
without applied magnetic field corresponds to the ver- 
tical axis. Then the two regimes R or RK are possible 
depending whether rflje > O -1 or not. Without both 
rotation and applied magnetic field a K type of turbu- 
lence is found. 

From our analysis we note that inertial regimes of type 
AK or RAK are never possible. On the other hand in- 
ertial regimes of type KA, A, or K are possible provided 
the forcing scale is sufficiently small. In Fig. [1] it corre- 
sponds to begin the spectra at a larger wave number. For 
Pm < 1 the inertial range of the kinetic energy spectrum 
prolongates at scales smaller than k^ with either an R, 
K or RK spectrum. 

III. SHELL MODEL 
A. The model 

The equations of MHD turbulence for an incompress- 
ible fluid embedded in an external uniform magnetic field 
Bo and subject to rotation H write 

du/dt + (u- V)u- ((B + Vj4 ) • V)B + 2£2 x u 

= ^V 2 u + F-VP f (7) 

dB/dt + (u-V)B-((B + v A )-V)u = r ? V 2 B (8) 

V • u = V • B = (9) 

in which va = Bo / J~\xp is the Alfven velocity (where 
fi and p are respectively the fluid magnetic permeability 
and density) and B is given in unit of V A = \v A \- The 
total pressure Pt = P + 6 2 /2 is a functional of u and B 
owing to the incompressibility condition (|9]). The forcing 
F insures the fluid motion. 

From (J7J (HJ) ([9]) we derive the following shell model 

U n = ik n [Q n (U,U)-Q n (B,B)} 

+ ik n V A {t)B n + %Sl(t)U n - vklU n + F n (t), (10) 

B n - ik n [Q n (U,B)-Q n (B,U)} 

+ ik n V A {t)U n - nk 2 n B n , (11) 
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FIG. 1: Possible inertial regimes of energy spectral density in rotating MHD turbulence for Pm = 1. The capital letters R, A 
and K denote a rotating, Alfven or Kolmogorov turbulence. 



Q= v/ Irj 



{sir,) 




has the advantage that helicities are much better defined, 
like those based on helical wave decomposition (20l - [2^ |. 
It has been introduced in its hydrodynamic form to study 
spectral properties of helical turbulence [23j], and in its 
MHD form to study cross- helicity effect on cascades [24| ■ 
The parameter A is the geometrical factor from which 
the wave number is defined k n = fcoA". As explained 
in [l3j an optimum shell spacing is the golden number 
A = (1 + \/E)/2. The terms involving fl and Va were al- 
ready introduced in several previous papers dealing with 
either rotation (2^, [2(| or applied magnetic field [23, [III . 



FIG. 2: The four possible turbulent inertial regimes given in 
the map (Va, 



where 

Qn(X, Y) = A 2 (AT„ + iy, l+ i + X* +1 Y* +1 ) - X r n _]Y n 
—X n Y n _ 1 + i\(2X n Y n _ 1 + X n+1 Y n+1 — X n+1 Y n+1 ) 
+Xn-iXn-i + X*_ 1 Y*_ 1 - X 2 (X 1 n+l Y n + X n Y t r l+1 ) 

+iA(2AT*i^ +1 + X^ l _ 1 Y J l_ 1 - X^ l _ 1 Y^_ 1 )[12) 

represents the non linear transfer rates and F n the tur- 
bulence forcing. This model is based on wavelet decom- 
postion Compared to other shell models [l7l - [l9| it 



B. Conservative quantities 



Expressions for the kinetic energy and helicity, Ejj and 
Hjj, magnetic energy and helicity, Eb and Hb, and cross 
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TABLE I: In 3D MHD turbulence, conservation properties of 
total energy E, cross-helicity He and magnetic helicity Hm 
depending on global rotation Q and applied field Va- 



helicity He, are given by 



Ejj 
Hu 
E B 
H b 



U, 



2 

n | j 



(13) 



J2 E u{n), Eu{n) 

n 

J2 H u(n), Hu{n) 

n 

Y, E B(n), E B (n) = ±\B n \ 2 , 

n 

£ff B (n), H B (n) = l -K\{B* n f - B 2 n ),(16) 



M( U n?-Ull (14) 



(15) 



J2Hc(n), H c (n) = -(U n B* + B n U*). (17) 



In the inviscid and non-resistive limit (y = 77 = 0), 
the total energy E — Ejj + Eb, magnetic helicity and 
cross helicity must be conserved (E = Hb = He = 0). 
Here with the additional Coriolis and Alfvenic terms the 
properties of conservation are not necessarily satisfied. 
A summary of theses properties is given in table U for 
3D MHD turbulence. In the case of pure hydrodynamic 
turbulence (without magnetic field) the kinetic energy 
and helicity must be conserved {Ejj = Hjj = 0) even 
with Coriolis forces. 



D. No injection of cross-helicity 

In addition it is important to control the injection of 
cross-helicity as was shown in (24[. Indeed any spuri- 
ous injection of cross-helicity may lead to a supercorre- 
lation state where U n « B n implying equality not only 
in intensity (as in equipartition) but also in phase. In 
that case the flux of kinetic energy is depleted, imply- 
ing an accumulation of energy at large scale and steeper 
spectral slopes. In order to compare the results to the 
phenomenological approach wc impose the injection of 
cross-helicity to be zero. For that we could use the forc- 
ing 



\F„ 



±i-. 



B,, 



\B n 



(18) 



where again the sign is randomly changed after each pe- 
riod of time tF- This forcing is however ill-defined as 



1. To fix this problem we use 



the following forcing 



F ni 



a + C 



(19) 



with £ 



\Bl F \ 



in which ip is a phase randomly 
changed after each period of time tp, and a an additional 
parameter. In the case ( ^> a, (fT8|) is recovered, and the 
phase of F np is mainly determined by the phase of B np 
so that it corresponds to zero injection of cross-helicity. 
In the case (<a the phase of F nF is controlled by the 
random phase ip. Since B riF is small there is no cross- 
helicity injection too. The value a = 10~ 6 provides a 
robust forcing with always a low level of cross-helicity. 



E. Dissipations 



C. Time-scales 

In (fTU)) and pi[) the forcing F nF (t) (applied at some 
scale k~ F ), the global rotation f2(t) and the applied field 
VA{t) have constant intensities |F nF |,fi and Va- Only 
their sign may change after a period of time tp, in and 
ty A , the probability of changing from one period to the 
next being random. Such a trick allows to control the 
two characteristic times sa tn and ry A sa ty A . In the 
simulations we take tn — 1/51 and ty A = l/(k n „ Va)- It 
is in same spirit than the one used in [28| and [25| though 
much simpler. Incidentally the random change of sign of 
0(t) insures that there is no injection of kinetic helicity 
on average. Taking a random sign in F nF (t) we insure 
that the forcing intensity satisfies \F nF \ w \f2ejtp. It 
is also important that tp is the shortest among all other 
characteristic times of the problem Tk,tq and ry A (and 
of course r n ). We choose tp < min {tk, tq, ry A }. 



We define the dissipation of U and B at scale k n by 
Du(k n ) = vkl\U n \ 2 and D B {k n ) = rik 2 n \B n \ 2 . From the 
phenomenological formalism above we expect the total 
dissipation to be equal to the injection rate of energy at 
the forcing scale e v + e n — e, with e„ = J2 n E > u{k n ) and 
c, ; = Yin E>B{kn)- Equivalently in pure HD we would 
have e v = s. On the other hand the ratio of both dissi- 
pations p = e v /e u cannot be predicted. It can only be 
calculated numerically. 



IV. RESULTS 
A. Spectra for Pm = 1 

In Fig. [3] the spectra are plotted for v = 10~ 7 and 
Pm = 1 in the three cases Va = 0, fl = 0, and VaQ 7^ 0. 
For Va = 0, the horizontal and k 1 ^ 3 dashed lines disclose 
a RK regime. For il = 0, the dashed line k 1 ^ 6 disclose 
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100 




a KA regime. For VaQ 7^ 0, the A: -1 / 2 and horizon- 
tal dashed lines disclose a RA regime. In each case the 
transition between two power laws is rather smooth and 
occurs over a scales range of about two orders of magni- 
tude. 

For £ m 1 and taking the numerical values for fi, Va 
and 77 given in Fig. [2] we find that the three sets of spec- 
tra found with the shell model belong indeed to the three 
parts RK, (R)KA and RA of Fig. O We tried to track the 
transition from one part to the other, varying and Va- 
It is however not possible to handle it numerically as the 
spectral slopes are not so well defined at the neighbor- 
hood of the frontiers delimiting the four parts of Fig. [5J 

B. Spectra for Pm < 1 
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FIG. 3: (Color online) Normalized spectra for v — 10 7 
and Pm = 1. Curves (a) are shown for Va — and Q = 
12.5,25,50,10,200,400,800,1600 (from right to left, from 
darker to lighter). Curves (b) are shown for Q = and V A = 
0.16, 0.32, 0.64, 1.28, 2.56, 5.12, 10.24, 20.48 (from left to right, 
from darker to lighter). Curves (c) are shown for (V A ,Sl) — 
(0.16, 800); (0.32, 400); (0.64, 200); (1.28, 100); (2.56, 50); 
(5.12, 25); (10.24, 12.5); (20.48, 6.25) (from left to right, from 
darker to lighter). 



In Fig. 0] the kinetic and magnetic spectra are plotted 
for v = 1CP 7 and several values of Pm, for the three 
previous cases. 

For Va = (a,b) increasing Pm decreases the mag- 
netic dissipation scale while the viscous scale is not sig- 
nificantly changed. This is in agreement with a simple 
Kolmogorov phenomenology [JJJ , the ratio of dissipation 
scales being given by k v /kr, oc Pm -3 ' 4 . For Pm > 10 -2 
the effect of rotation is visible in the spectra flatness. At 
smaller values of Pm it is however difficult to determine 
any slope at all. 

For f2 = and Va = 1-28 (c,d) both kinetic and mag- 
netic spectra are almost the same whatever the value of 
Pm. The effect of an applied magnetic field is to corre- 
late both fields as expected in Alfven waves. In particular 
the dissipation scale is governed by the magnetic diffu- 
sivity, with k v ~ k n . The same conclusions arc found 
for (V A ,fl) = (0.32,400) (e) and (V A ,Q) = (20.48,6.25) 
(f). In these two cases the horizontal slopes are due to 
rotation (e) and applied magnetic field (f). 

We note that for ft = and V A = 1.28 (d) the nor- 
malized curves are not horizontal. They correspond to 
spectral energy density slopes between /c~ 5 / 3 and fc~ 3 / 2 . 
The latter is obtained for values of V A about ten times 
larger. 

C. Dissipation ratio 

In Fig. [5] the ratio p = e n /e v is plotted versus Pm for 
V A = H = (a), V A = (b), Q = (c) and V A Q ^ 
(d). In the limit Pm — > the dynamo action does not 
occur, implying p — > 0. For Pm = 1 both kinetic and 
magnetic spectra are identical, implying e v — e v = e/2, 
and then p = 1 . We always find an intermediate value of 
Pm for which p reaches a maximum. This is related to a 
super-equipartition state in which the magnetic energy is 
higher than the kinetic energy at large scales. Varying Va 
and SI we find that this maximum value can increase by 
several orders of magnitude and that it does not occur at 
the same Pm. For the two last cases an asymptotic curve 
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FIG. 4: (Color online) Normalized kinetic (full) and magnetic (dashed) spectra for v — 10 7 : (a) (Va,£1) = (0,100), (b) 
(Va.O) = (0,1600). (c) (Va,0) = (0.08,0), (d) (V A ,ty = (1.28,0). (e) (V A ,0) = (0.32,400), (f) {V A ,Sl) = (20.48,6.25). For 
each set of curves for Pm = 10~ 5 , 10 -4 , 10 -3 , 10 -2 , 10 _1 , 1 (from lighter to darker). Note, that the kinetic and magnetic spectra 
are superposed in the case (f). 

p = 0(Pm _1 ) is obtained for large values of V A . This is V. DISCUSSION 

a direct consequence of the equipartition regime | U n \ m 
\B n \ obtained at any scale (see Fig.|4|). In that cases the 

definition of p directly implies the scaling 0(Pm _1 ). For Pm = 1 both approaches, phenomenological and 

shell model, give consistent results in terms of inertia 
regimes. They are controlled by the shortest time-scale 
corresponding either to rotation, applied magnetic field, 
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FIG. 5: Dissipation ratio versus Pm. On panel (a) Va = = and the full curves from right to left correspond to v — 10 5 , 
1CT 6 , 1(T 7 , 10 -8 . The dashed curves from bottom to top correspond to i) = 1/4, 1/8, 1/16, 1/32, 1/64, 1/128, 1/256. On panel 

(b) Va = 0, v = 1(T 7 and the curves correspond to = (full curve), 100 (dashed), 400 (dotted), 1600 (dot-dashed). On panel 

(c) f2 = 0, v = 10~ 7 and the curves correspond to Va = (full curve), 0.08 (dashed), 0.32 (dotted), 1.28 (dot-dashed). On 
panel (d) v = 10" 7 and the curves from bottom to top correspond to [Va,Q) = (0,0), (0.32, 400), (2.56,50) and (20.48,6.25). 



inertia, or a combination of them. For Pm < 1 the mag- 
netic dissipation occurs at a scale larger than the viscous 
scale implying that the different regimes are not so easy 
to discriminate. However for a sufficiently strong applied 
magnetic field both kinetic and magnetic energy spec- 
tra are merged, implying a strong increase of the viscous 
dissipation scale. Whether this is due to our isotropic 
assumption is not clear and cannot be answered with 
our models. A consequence is that, for a strong ap- 
plied field, the ratio of magnetic to kinetic dissipation 
scales like O(Pnr') and can reach very high values for 
Pm <C 1. Without applied field, this ratio is also maxi- 
mum for some value of Pm -C 1 , depending on the fluid 
viscosity and global rotation. 
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